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Abstract 

CD . There now exists in the hterature two different expressions for the phase shift of a matter- 
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CN '_ wave interferometer caused by the passage of a gravitation wave. The first, a commonly accepted 

^ ■ expression that was first derived in the 1970s, is based on the traditional geodesic equation of 
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motion (EOM) for a test particle. The second, a more recently derived expression, is based on the 
geodesic deviation EOM. The power-law dependence on the frequency of the gravitational wave 



^ ■ for both expressions for the phase shift is different, which indicates fundamental differences in the 

On . physics on which these calculations are based. Here we compare the two approaches by presenting 

^S. '. 8- series of side-by-side calculations of the phase shift for one specific matter-wave-interferometer 



configuration that uses atoms as the interfering particle. By looking at the low-frequency limit of 



O . the different expressions for the phase shift obtained, we find that the phase shift calculated via 

^H . the geodesic deviation EOM is correct, and the ones calculated via the geodesic EOM are not. 
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I. INTRODUCTION 

In two recent papers [l|, |2| , a new calculation of the phase shift of a matter- wave inter- 
ferometer caused by the passage of a gravitational wave was presented. This calculation 
was based not on the usual geodesic equation of motion (EOM), but rather on the geodesic 
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'orm that is very much 
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deviation EOM. Contrary to expectations, these phase shifts 
different from those calculated previously in the literature 
direct comparison of the various expressions for the phase shift is difficult because of the 
different configurations of the matter-wave interferometers used for their calculations, it was 
noted in |1] that the power-law dependence on the frequency of the gravitational wave, in 
both the low- and high-frequency limit, of the expressions calculated in 111 yl are at times 
d..e**any oppos.e .o .Kose calculated .0^03030. C„.c,l%, we expect 
that the differences between these expressions are not due simply to differences in the con- 
figurations of the interferometers chosen, but are rather caused by fundamental differences 
in the physics on which their derivations are based. 

In this paper our objective is to demonstrate that these fundamental differences do exist 
by presenting a series of calculations that follows the two basic approaches in the literature 
used to determine the phase shift caused by a gravitational wave. This will be done for 
the same configuration of matter-wave interferometer, and will allow us to compare the 
different expressions we obtain for the phase shift on an equal footing. We then determine 
their validity by comparing each of them, in the appropriate limit, with the well-known 
expression for the phase shift caused by stationary sources of spacetime curvature |5l.llCilll|. 
The focus in this paper will be on a theoretical study of the underlying physics: we leave 
issues of the practicality and feasibility of constructing the interferometer to [l|,y]. 

The various approaches taken to determine the phase shift of a matter- wave interferometer 
caused by a gravitational wave can be divided into two distinctly separate classes. The class 

transverse-traceless (TT) coordinates (see Sec. 9.2.3 of [1^); either explicitly or implicitly. 
The TT coordinates are one specific choice of reference frame where the gravitation wave 
is in the TT gauge. In this frame the number of components of the deviations of the flat 
spacetime metric h^^ caused by the gravitational wave equals the number of physical degrees 
of freedom of the gravitational wave. Consequently, d'^h^i,, /iq^ = 0, and h'^ = 0. In addition, 



the mirrors and the beam sphtters are at fixed coordinate values in these coordinates, and 
they do not move. These coordinates have been used to study the properties of hght- 
wave-interferometry-based detectors of gravitational waves as well |13]. To calculate the 
phase shift for matter-wave interferometers, it is combined with either the geodesic EOM 
(and ,t. action). ,f the calculat.on . done at the ,uantun,-,neehaniea. .eve. fl Q. o. t.. 
Lagrangian (and the WKB or stationary phase approximation) for a quantum field in curved 
spacetime, if the calculation is done at the quantum field-theoretic level p, ly, LD, |8| . Not 
surprisingly, irrespective of the level of sophistication of the approach used, the general 
expressions calculated for the phase shift have the same final form. 

The other class of approaches to calculating the phase shift is taken by [l|, ^L and is done 
in Thome's "proper reference frame". They follow the classical analysis in |l2l. Il4l. Il5| for 
light-wave-based interferometers. In this frame, the positions and velocities of test particles 
are measured relative to the worldline of an observer, and follow the motion of the observer in 
.paceti^e. ,n .ach the .a.e wa. that the Fe..,-ao.n,a. U . Fe.n„-Wa..e. coo.d.nate. Q. 
and the recently constructed general laboratory frame 1^ do. Consequently, the motion 
of the test particle in this frame is not described by the usual geodesic equation of motion 



(EOM), but rather by the geodesic deviation EOM (see Eq. (2) of 1^)) which was first 



applied to the analysis of physical systems in Il9l l2Ci . As in |l^], the TT gauge was taken 
for the gravitational wave. The action (see [l8|, |2l|, |22|) for the geodesic deviation EOM 
was then used in conjunction with the quasi-classical approximation of the Feynman-path- 
integral formulation of quantum mechanics. This approximation is the Lagrangian form of 
he WKB approximation in the Schrodinger representation of quantum mechanics, and in 
l|,l2| it was applied to the calculation of the phase shift of the two matter- wave-interferometer 
configurations. 

The phase shift of a matter-wave interferometer is a physically measurable quantity. 
From the principle of general covariance, it should not matter whether the TT coordinates, 
or whether the proper reference frame is used in its calculation |l2L l23l |. Indeed, Thorne 
states the following: 

"If L « X then the [gravitational wave] detector can be contained entirely 
in the proper reference frame of its center, and the analysis can be performed 
using non-relativistic concepts augmented by the quadrupolar gravity-wave force 
field (3), (5). If one prefers, of course, one instead can analyze the detector in 



TT coordinates using general relativistic concepts and the spacetime metric (8). 
The two analyses are guaranteed to give the same predictions for the detector's 
performance, unless errors are made [emphasis added]." (From page 400 of [12] •) 

Evidently, errors have indeed been made, as we shall see, since the resulting calculations 
often contradict each other for the same interferometer. The power-law dependence on 



frequency of the phase shifts calculated in 
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using TT coordinates are 



fundamentally different from that of the phase shifts calculated in |l|, |2] using the proper 
reference frame. Since theDrinciple of general covariance must hold^this means that either 
(a) the calculate,. . Q^ a.e wong. (b) the calculat.o,. .. BflfiBflaa are w.ong. 

or (c) all of the calculations of the phase shifts are wrong. 

To determine which one of these three possibilities is correct, we will present side by side 
calculations of the phase shift caused by a gravitational wave for one specific matter-wave- 
interferometer configuration using the following three approaches. The first calculation uses 
the quasi-classical approximation in the TT coordinates and geodesic action, and is done 
at the nonrelativistic, quantum-mechanics level. It follows, as much as possible, the line of 
analysis given in |Jj. The second calculation is also done at the nonrelativistic, quantum- 
mechanics level using the same approximation, but now in the proper reference frame, and 
uses the geodesic deviation action. It follows the line of analysis given in \^ . 

These two calculations would seem to exhaust the two different classes of approaches 



currently found in the 
because of an error in 



iterature. That a third calculation of the phase shift is needed is 
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term in the action for the atom in 



]. Either because changes in the kinetic energy 
were neglected, or because the turning points of t 



3omts oi the 



atom's phase at the interferometer's mirrors in the WKB-type of approach in 
were omitted, the end result is that important boundary terms are not included in the final 
expre..o„ fo. the phase shift in M fl B fl B 3. Fo. the saUe of eo,np.ete.e., and 
to compare the calculations in [l|, |2 1 with these previous calculations in the literature, this 
necessitates a third calculation of the phase shift, which will result in yet another expression 
for the phase shift based on the geodesic EOM. 

We find that all three expressions for the phase shift calculated here are different from 
one another. Although there are currently no direct experimental checks that can used to 
determine which one of these three expressions for the phase shift is correct, an indirect 
experimental check can be used by taking the low-frequency limit. 



All three expressions for the phase shift depend on the transit time of the atom through 
the interferometer, and the period of the gravitational wave. If the transit time of the atom 
is much shorter than the period, then, with respect to the atom, the gravitational wave 
does not oscillate appreciably as it traverses the interferometer. The gravitational wave will 
contribute a local Riemannian curvature to the spacetime that is effectively static. The phase 
shift of matter-interferometers caused by static (as well as stationary) sources of curvature 



is well established 
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ll|, and has been shown to be proportional to components of the 



Riemann curvature tensor. Just as importantly, this phase shift can be calculated using only 
Newtonian gravity [Eq. (2.6) of (lO]]. The phase shift due to the Newtonian gravitational 
potential — through the acceleration due to Earth's gravity g — has been measured, first by 
CoUela, Overhauser and Werner [2^, and more recently by Chu and Kasevich |25|. These 
more recent measurements of the phase shift induced by g are sensitive enough that effect 
of spatial variations in g — which are proportional to the static Riemann curvature of the 
spacetime caused by the Earth — on the phase shift can be discerned. Gradients in g were 



seen throug 
was varied 



1 changes in the phase shift of an atomic fountain when its vertical position 



2y|. More recently, they were seen through a differential measurement of the 
phase shifts of two atomic fountain matter- wave interferometers operating in tandem J27|. 
Therefore, while the direct measurement of the phase shift of a matter-wave interferometer 
caused by a static Riemann curvature tensor has not been made, its discovery is expected. 



Indeed, indications are that t 
Earth has already been seen 



the phase shift induced by the static Riemann curvature of the 

As a consequence of the above, in the limit of low-frequency gravitational waves the 
phase shift calculated for the matter-wave interferometer must approach the well-known 
static result, and be proportional to the curvature tensor. We will see that in this limit, 
only the phase shift calculated along the lines of [l|, |2] has the correct form. The other two 
do not, and must be wrong, since they do not yield the correct static limit. 

Our focus here is on the phase shift that a gravitational wave will induce on a matter-wave 
interferometer. We will thus neglect the effects of all other gravitational effects — such as the 
acceleration due to Earth's gravity, and the local curvature of spacetime due to stationary 
sources such as the Earth, Sun and Moon — in addition to the phase shift due to the Sagnac 
effect caused by the Earth's rotation. We have shown in [l8| that for these types of stationary 
gravitational effects the description of the motion of test particles in the general laboratory 



frame — of which the proper reference frame is a special case — reduces to the usual geodesic 
EOM description, and there is no controversy. The phase shifts caused by these stationary 
effects calculated in either the proper reference frame or the TT coordinates will be the 
same. 

We will in this paper use the same notation the is commonly found in the literature 
111, y, y, y, y, y, li y, bl- in particular, it follows the notation in |jj. While this has the 
benefit of consistency and simplicity, it means that we do not separate with notation the 
coordinates used for the TT coordinates from the coordinate used for the proper reference 
frame. The coordinates used in the two reference frames have different physical meanings, 
which we will make clear in the next section. These differences are outlined in detail in [l8[ . 

II. PHASE-SHIFT CALCULATIONS 

In this section we will present two of the three calculations of the phase shift mentioned 
in the introduction. The first calculation will be done in the TT coordinates, and will be 
based on the geodesic EOM. The second calculation will be done in the proper reference 
frame, and will be based on the geodesic deviation EOM. A third calculation, following the 
approach in the literature, will be done in the next section where comparisons between the 
different results of the calculations are made. All three calculations are done for gravita- 
tional waves in the long-wavelength limit for the matter-wave interferometer configuration 
shown in Fig. ^ this configuration is described in detail in appendix A. To calculate the 
phase shifts, we will use the quasi-classical approximation of the Feynman path integral 
representation of quantum mechanics. This approximation is the Lagrangian version of the 
WKB approximation, or stationary phase approximation, in the Schrodinger representation 
of quantum mechanics, and will be reviewed in appendix A as well. The presentation in 
this section is done in detail to ensure that all the underlying assumptions, approximations, 
and subtleties in our analysis are readily apparent. 



A. Geodesic-EOM-based calculation of the phase shift 

1. Classical Dynamics in TT Coordinates 

Consider a gravitational wave in the long- wavelength limit with an amplitude hfj,i^{t) 
incident perpendicularly to the plane of the interferometer in Fig. ^ As usual, Greek indices 
run from to 4, and Latin indices run from 1 to 3. The gravitational wave is represented by 
the tensor h^iy as a perturbation on the flat spacetime metric rj^iy, so that the total metric 
of the spacetime is g^^ = tj^^ + h^^,. We have chosen the signature of g^y to be (—1, 1, 1, 1). 
In the TT-gauge, a coordinate system is chosen so that in these coordinates, d'^h^^, = 0, 
Hq^ = 0, and /i^ = 0; the only nonvanishing components of h^^ are thus hij. 

One choice of coordinates where the TT gauge is realized is the TT coordinates (see [12| 
and the general analysis in [29|). In these coordinates the action for the atom with mass m 
used in Fig. ^is the usual geodesic action in the nonrelativistic limit, 

( — dflU^'''' I ^^ ^ -mcH + 2"^ / ["''"^i^ij + ^ij)] dt, (1) 

where the subscript "geo" stands for "geodesic". Since the rest mass mc^ only shifts the 
Lagrangian for the atom by a constant, we drop this term from Sgeo- The nonrelativistic 
Hamiltonian for the motion is then, 

which results in the usual geodesic equation of motion 

dt^ ~ ^^ '' dt ~ "^ dt • ^^^ 

Here F" is Christoffel symbol for the spacetime; for gravitational waves in the TT gauge, 

To calculate the phase shift of Fig. ^ we follow the description in appendix A, and 
consider the phase shift due to the rth atom emitted from the atomic source, which at 
a time tj. is diffracted by the initial beam splitter. There are two possible classical paths 
for the atom through the interferometer. In the absence of gravitational waves, these are 
straight-line paths xo{t) (see appendix A), which will be shifted when a gravitational wave 
is present. These shifts are expected to be very small, however, and we can use perturbation 
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FIG. 1: Schematic of the matter-wave interferometer configuration considered in this paper with 
a gravitational wave incident normal to the plane of the interferometer shown. Only the x polar- 
ization contributes to the phase shift; the + polarization does not. Diffraction orders other than 
n = ±1 are not used, and are left out for clarity. 

theory to solve for them. We thus take x{t) = Xo(t) + Xi(t), and v(t) = Vo(t) + Vi{t), where 
xi(t) and vi(t) are perturbations to the atom's path xo{t) and velocity vo(t), respectively; 
we will keep only terms linear in vi{t). To first order in hij, Eq. becomes 



dt ~ ^^'''■ 



(4) 



Once Xo{t) is determined, fi(t) is solved for by dividing the trajectory of an atom into two 
parts: from the initial beam splitter to the mirror, and then from the mirror to the final 
beam splitter. The two paths are then joined with the appropriate boundary conditions at 
the mirror. 

We begin by determining Xo(t). 



2. Trajectories in the absence of gravitational waves 

To determine the unperturbed path of the atoms through the interferometer — which is 
the same no matter which approach is used to calculate the phase shifts — consider first the 
path of the atom that is diffracted into the n = +1 diffraction order. We will assume that 
the velocity of the atoms is high enough that we can ignore the effect of Earth's gravity on 
the atom. It is then clear that for t^ < t < t^ + T/2 — during which it travels between the 
initial beam splitter and the first mirror in a time T/2 — the path of the atom is 

V0x=V\\, Xo{t) =Vl\{t-tr), 

VOy = V±, yo(t) = f±(t - tr), (5) 

where v\\ and v± are x and y components of the velocity of the atom immediately after it 
leaves the initial beam splitter. At tr+T/2 the atom is reflected off the hard-wall mirror, and 
from the general expression of the boundary condition for the atom at the mirror Eq. ()A2|) . 
Vo,{tr + T+/2) = vo,{tr + T-/2) = v\\, while Voy{tr + T+/2) = -voy{tr + T-/2) = ~v^. (The 
superscripts ± denote lime-^o^r + T/2 ± e respectively.) For the subsequent time interval 
tr + T/2 < t < tr + T, the atom's path is 

T;0:r = f||, Xo{t) =Vll{t-tr), 

voy = -v±, yo{t) = L±- v^{t - tr), (6) 

during which it travels between the mirror and the final beam splitter. It is clear that L_j_ = 
Vj_T and Ly = f ||T, and T is the total transit time of the atom through the interferometer. 
The unperturbed path for the atom in the n = —1 diffraction order is similar to Eqs. (jSJ 
and (jBI), and is obtained by taking v± -^ —v± and L± -^ —L±. 

3. Boundary conditions 

To determine the boundary conditions for vi in Eq. (j3)), we note that from Eq. (jS)) the 
acceleration of a test particle for the geodesic EOM is proportional to its velocity. As such, a 
test particle which was initially at rest will stay at rest, even in the presence of a gravitational 
wave. This led to the following conclusion by Stodolsky: 

"In this discussion [of matter-wave interferometry] we have assumed that the 
parts of the apparatus [interferometer] stay at fixed coordinate values. This can 
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be achieved if these parts are free masses, since a metric of the gravity wave type 
(with no "0" components) does not affect the positions of nonmoving bodies." 
(From page 399 of [^].) 

Thus, in the TT coordinates, the positions of beam sphtters and the mirrors do not change 
when a gravitational wave passes through the system. 

The velocity of the atoms travelling between the supersonic source and the initial beam 
splitter, on the other hand, will be affected by the gravitational wave. Since the interfer- 
ometer in Fig. H is balanced, only differences in the atom's trajectory along the two paths 
through it will matter. We can thus neglect any shifts in the atom's velocity before it reaches 
the beam splitter. And, since the initial beam splitter is not moved by the gravitational 
wave, voxifr) = v\\ while voy(tr) = v±; thus we take viifr) = 0. 

For the boundary conditions at the mirror, we note that like the beam splitter, the 
positions of the mirrors in the interferometer also do not change when a gravitational wave 
is present. Thus the boundary conditions for Vi{t) across the hard- wall mirrors is 

vUtr + T+/2) = vUU + T-/2), viyiU + T+/2) = -v.yiU + T' /2). (7) 

Integrating Eq. (Q is now straightforward, and for tj. < t < t,. + T/2, 

Vlx{t) = -V\\K^it) - V±Ky{t) + Vllh^xiU) + V±h,^y{tr), 

^lyit) = -V\\hxy(t) - V_Lhyy(t) + V\\h^y(tr) + V±hyy{tr), (S) 

while for U + T/2 <t<tr + T, 

Vlx{t) = -V\\hxx{t) + V±hxy{t) + V\\hxx{tr) " 2v±hxy{tr + T/2) + V±hxy{tr), 

Vly{t) = -V\\hxy(f) + V±hyy(t) - V\\h^y{tr) + 2v\\hxy{tr + T/2) - V±hyy{tr). (9) 

4- The phase shift 

We now calculate the action for the rth atom as it travels through along the upper path 
through the interferometer. To first order in hij, 

Sfjx{tr),x{tr + T);^^,] = ^mj^ {[voit)+v^it)f + vUh^J}dt, 

-ml Vq dt + m <vo- vi{t) + -fX^ij \ dt,{10) 
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where the superscript "cl" stands for "classical" , and denotes that the integration is over a 
classical path through the interferometer. The Vq term in Eq. (fTUI) is the same for both paths 
through the interferometer, and will not contribute to the phase shift. It can be neglected. 
The last term VQvlhij can be trivially integrated using Eqs. (0) and (jH)), while the Vq ■ Vi term 
can be integrated using Eqs. (jSj) and Q- We then get 



S^J°[x{tr),x{tr + T); 7+1] =T [vJK^{tr) + vlhyy{tr)] 

1 

~ 2 



tr+T r 



vjh^xit) +v\hyy{t) 



dt 



+ 2v±v\\T Ky{U) - Ky{tr + T/2) 



•'tr+T/2 



Vl_V\\ 



hxyit) dt 



tr+T 



hxy{t) dt 



:ii) 



'tr J tr+T/2 

where the tilde denotes the fact that we have dropped the v"^ term from the action. We 
get the action Sg^J[x{tr) ^ x{tr + T);7_i] for the atom if it took the lower path through the 
interferometer by taking v±_ -^ —v± in Sg[^[x{tr),x{tr + T);7+i]. The first two terms in 



Eq. (fTT|) do not change, while the last two terms flip sign. Thus, from Eq. ()A1|) . the phase 
shift of the atom is 



2/72 



h 



itxyyj'r ' -'- / ^ ) i^xyy^r) 

T 



•'tr+T/2 



+ 



2^2 



hxy{t) dt 



tr+T 



hxyit) dt 



tr+T/2 



After taking the Fourier transform, 

duj . 



hxy{t) = 



27r 



hxi^)e 



-iujt 



geo 



it) 



duo 

2^ 



A0geo(u;)e' 



-iujt 



(12) 



(13) 



the calculated phase shift for Fig. ^ using TT coordinates simplifies to 



m 

2h' 



A(f)gU^) = - ^L\\L^iuoK{uj)e-'^^'^Fg,o{ujT), 



geo 



uT) 



sine (u;T/4) 



4e^'^^/^sinc(u;T/4), 



(14) 



where sine (x) = sin(x)/x. The first result of our comparison of phase-shift calculations for 
the interferometer shown in Fig. ^is given in Eq. ^^. The derivation of these equations is 
based on the geodesic equation of motion using the TT coordinates. We will find that the 
phase shift calculated in all three approaches in this paper will have this same overall form, 
and will differ only in the resonance function Fg^o- 
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B. Geodesic-deviation-EOM-based calculation of the phase shift 

Reference frames centered on the worldline of an observer, such as the proper reference 
frame and general laboratory frame, are less commonly used than the TT coordinates used 
in the previous subsection. We therefore begin the calculation of the phase shift for Fig. Q 
in the proper reference frame with a brief review of this frame. 

1. The proper reference frame 

How the proper reference frame is constructed is described as follows: 

"Consider an observer (freely falling or accelerated, it doesn't matter so long 
as the acceleration is slowly varying). Let the observer carry with herself a 
small Cartesian latticework of measuring rods and synchronized clocks (a 'proper 
reference frame' in the sense of Sec. 13.6 of MTW, with spatial coordinates x^ 
that measure proper distance along orthogonal axes)." (From page 339 of ^12].) 

In this frame the coordinate a;*(t) is the position of a nonrelativistic test particle travel- 
ling along its worldline relative to the worldline of the observer at her propertime t. The 
magnitude of x*(t) is then the distance between the worldline of the test particle, and the 
worldline of the observer at the same instant of time t. This frame is only valid when the 
test particle is moving nonrelativistically, and the size of the apparatus is small compared 
to the wavelength of the gravitational wave. The velocity f*(t) = x* is then the velocity of 
the test particle measured relative to the observer. 

In the absence of all other forces, both the test particle and the observer travel along 
geodesies in the spacetime, and thus their motion are each, separately, governed by geodesic 
EOMs. By subtracting the geodesic EOM of the observer from that of the test particle, and 
then expanding the resultant to first order in x*(t), the motion of the test particle as seen 

by the observer is governed by 

d x^ 

-^ = -RqjqX^. (15) 

This is the geodesic deviation EOM, and it holds as long as the local Riemann curvature 
tensor of the spacetime Rq,:q varies so slowly with x* that it can be approximated as a 
constant. This condition holds in the long- wavelength limit of the gravitational wave. 
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As observed in |12] for LIGO (the Laser Interferometer Gravitational-wave Observatory), 
if the characteristic size of the interferometer L is much smaller than the reduced wavelength 
A of the gravitational wave the interferometer is contained entirely within the proper refer- 
ence frame. One does not need to use TT coordinates. Like LIGO, the characteristic size 
of the interferometers considered in 111 ^ is also smaller than the reduced wavelength of the 

riQ 

gravitational wave [H, I2|, and thus the use of the proper reference frame was also valid. We 
shall assume that this is true for the interferometer shown in Fig. IT] as well. 

we wUI also follow Tho.e Q. and take B,, ^ -1^,2. This .elation between B,,, and 
/i* holds only when the "TT gauge" is taken. Because of the naturalness of its construction, 
it should be expected that the TT gauge — where the number of independent components 
of hij equals the number of degrees of freedom for the gravitational wave — can be taken for 
the proper reference frame as well as for the TT coordinates, but this has not been shown 
explicitly in the literature J30|. As we shall see in appendix B, this can be done explicitly 
in the general laboratory frame. It can thus be done in the proper reference frame since it 
was shown in [l^ that the general laboratory frame reduces to the proper reference frame 
in the long-wavelength limit for the gravitational waves. 

2. Classical Dynamics in the proper reference frame 

For our calculation of the phase shift in the proper reference frame using the geodesic 
deviation EOM, the observer is fixed on the initial beam splitter, and the motion of all 
objects, and thus the phase shift, is measured with respect to it (see Fig. d). In this 

n n 

reference frame, the action for a test particle is (see [1^ or |22| for a derivation) 

Sgd = m ( ^^^ - ViX^'^oj j dt = -m [v^ - v'x^hijj dt, (16) 

where the subscript "GD" stands for "geodesic deviation", and we have dropped the mc^ 



term. As in Eq. (4) of [12] , we have taken the gravitational wave to be in the TT gauge in 
Eq. ()16|) . The resulting Hamiltonian (see also '21|) in this case is 

to lowest order in hij. It is straightforward to see from Eq. (fT?)|l that the EOM for the test 
particle is the geodesic deviation equation Eq. (jl5|) . We once again solve perturbatively for 
the classical trajectories, which are determined now by Eq. (fT3|). and take x{t) = Xo{t)+Xi{t). 
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3. Boundary conditions and mirror effects 

Unlike the geodesic EOM [Eq. Q], in the geodesic deviation EOM [Eq. (fT3|) ] the accel- 
eration depends on the position of the test particle, and not on its velocity. Consequently, 
when a gravitational wave passes through the system, it shifts the position of all parts of the 
interferometer except the initial beam splitter, where x* = due to our choice of origin. We 
can thus take Vi{t^) = as an initial condition for vi{t) once again. In addition, from the 
force lines drawn in Fig. ^ we see by symmetry that it is the x polarization that contributes 
to the phase shift, and not the + polarization. For this polarization, the motion induced on 
the final beam splitter by the gravitational wave also does not affect the overall phase shift. 

While the motion of the final beam splitter does not affect the phase shift of the atoms, 
the motion of the mirrors in response to the gravitational wave will. This motion changes 
the boundary condition for the atom at the mirrors from the simple expression Eq. ()A2|) . 

Let X denote the position of the mirror placed in the upper path of the interferometer, 
and let V denote its velocity. In absence of the gravitational wave the mirrors do not move, 
and we can take X = L\\/2 + Xi and Y = L±/2 + Yi, where Xi are fluctuations in the 
positions of the mirror due to the gravitational wave. Because of these fluctuations, the 
boundary condition Eq. (|7j) for Viy(t) at the mirror now becomes 

ViyiU + T+/2) = -v^y{tr + T+/2) + Viy{tr + T/2), (18) 

where Vy = Y = Yi. 

To determine the response Viy{t) of the mirror to the gravitational wave, we follow [l|, l^l, 
and model the mirror, and its connection to the frame of the interferometer, as a spring 
with a quality factor Q and resonance frequency ujq (see also section 37.6 of [15]). This 
frequency depends on the size of the interferometer as well as its material properties. From 
the equivalence principle, the motion X{t) of the mirror, like that of the atom, is also a 
solution of Eq. (fT3j) . but because the mirrors are connected to damped springs. 



Q"^ ' "'J"^^^^ ^L\\h,,+ ^. 



Xi + —Xi + co'gXi(t) — -L||/i2:x + -L±hxy, 



Y, + ^Fl + Uj'^Y.it) = h^hyy + h^h^y. (l9) 
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4- The Phase Shift 



The action for the atom along the upper path through the interferometer in the proper 
reference frame is 



Sajj[x{tr),x{tr + T];-f+i) = -m 



1 
^ —m 



tr+T 



tr+T 



[voit) + MW - ^o^hij \ dt 



Vn dt + m 



tr+T 



2 



1 ,•; 



ClJjn 



vu-^xih,A^dt, (20) 



where once again "cl" stands for "classical" . The Vq term is once again the same for both 
paths through the interferometer, and can be neglected. Performing an integration by parts, 
and using Eq. (fT3j) . we get 



1 X 



-S^l''[x{tr),x{tr + r);7+i] =X^(t) Vu{t) - -xi{t)h,,{t) 
lit \ Zi 

1 . 



1 . 



^x\{t)[vu{t)-^x^^{t)K,{t) 



^ l'tr+T/2 _ -^ 

+ - I x^vlhij dt + - 



tr+T/2 

tr+T 

tr+T/2 
tr+T 



tr+T 



Of] a ^^* 



(21) 



We integrate by parts once again, and then use Eqs. (0), and (jH)) to get 



1 X 



-S,^^[f(t,),x(t, + T);7+i] =L|| <^ Vi^{tr + T) - -L\\h^^{tr + T) 
m I 



+ -Lx {v^y{tr + T-/2) - v^y{t, + T+/2)} 



'rlfM<^)h^m 



tr+T/2 



+ 7^xl{t)vi{t)h,,{t) 



tr+T 



tr+T/2 



tr+T 



v^vlhij dt. 



(22) 



Obtaining SQi^[x{tr)jX{tr + T);7_i] the usual way from SQj^[x{tr)jX{tr + T);7_|_i], we 



arrive at 



m 



,(+1)^ 



,(-1)^ 



A0GD(t.) = ^{L\\ nfiU + T) - v'f,'>{tr + T) - v^KyiU + T) 



+ v±Kyitr + T/2) 



+ L 



± 



'v[-^'\u + T-/2)+v[-'\tr + T-/2) 



--Lj 
2 ^ 



vi:'\u + T/2) + vl-'\u + T/2) 






ftr+T/2 ntr+T 

- 2f ||f_L I / h^yit) dt- j Ky{t) dt 

'tr Jtr+T/2 



(23) 
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where we have used the boundary condition Eq. (jTHj). The superscripts (±1) denote whether 
it is the velocity of the atom or mirror along the top [represented by the superscript (+1)] 
or the bottom [represented by the superscript (—1)] path through the interferometer. 
We note that to lowest order in hij Eq. (J15j) becomes 



dt ~ 2^^'''^^^- 



(24) 



We then take the component of Eq. ()24j) along the x direction for an atom travelling along 
the lower path, and subtract it from the same equation for the atom travelling along the 
upper path. Since Xq (t) = Xq (t), while i/q (t) = —i/q (t), we integrate and find 

v[t'\tr + T) - v[-'\tr +T)= V^{h,y{tr + T) - 2/l,.,(t,, + T/2) + h,y{tr)y (25) 

Similarly, if we take the sum of these equations along the y direction, we get 



,(+1), 



,(-1)^ 



1 



V[;'>{tr + T-/2) + v[/>{tr + T-/2) = -^\Ky{tr + T/2) - V\\\^Ky{tr + T/2) - KyiU) j. 



Consequently, 



2m I 1 



h 



AL 



Vi^'\U + T/2) + Vi~'\u + T/2) 



\Ky{U + T/2) + ^^^^^^'^^-^^^^'^^ 



1 



tr+T/2 



T 

tr+T 



hxy{t) dt- Ky{t) dt 

J tr+T/2 

Using Eq. (fT^ and the steady-state solution for Eq. (fT^ . then 



(26) 



virHt. . m . <-(. . m - -\l, I '^^-r^^^- 



-iu}{tr+T/2) 



(27) 



and we find that 



m 



A(Pgd{u;) = - —^L^iujh^{uo)e-'-^'^FGD{i^T), 



FcoiioT) =1 - 2e^'^^/Sinc (ujT/A) + 



sine {ujT/i) 



-1 2 1 2 

1 u'^ 



(28) 



2 cij^ — cUq + iuJUJo/Q 

The second result of our comparison of phase-shift calculations for the interferometer shown 
in Fig. ^is given in Eq. (j28p . The derivation of these equations is based on the geodesic 
deviation equation of motion using the proper reference frame. The phase has precisely the 
same form as Eq. (|TH): only the resonance function Fgd has changed. 
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Quantity 


Geodesic-equation based 


Geodesic-deviation-equation based 


Action 


mjdt[v'v^{6ij + hij)] /2 


m f dt 


v'^ - v^x^hij 


/2 


EOM 


p = -A}x^' 


x' = 'hi.xi/2 


Hamiltonian 


(p2 - hijpY) /2m 

_ _ 1 


p'^/2m + hijx'p^/2 

r " ^ 


Resonance Function 


sine (ujT/A) 


- 4e*^^/^sinc (coT/i) 


1 - 2e*'^^/^sinc {ujT/4) + 


sine (a;T/4) 








-J^ 1 {ij? - t^o + ^w^o/Q) /2 


Low-frequency Phase Shift 


-3mL|,Lxrg^(a;)/n 


-mL||L_LTi?ozOj/(w)/2?i 



TABLE I: A comparison between the geodesic EOM and the geodesic deviation EOM descriptions 
of test-particle dynamics, and the phase shifts calculated for Fig. fusing the two EOMs. 

III. THE LOW-FREQUENCY LIMIT, AND DETERMINING THE CORRECT 
PHASE SHIFT EXPRESSION 

A summary of the geodesic- and geodesic-deviation-based EOM approaches to the de- 
scription of the motion of test particles in the presence of a gravitational wave, along with 
the resonance functions for the phase shifts calculated with both approaches, is listed in 
Table HI It is readily apparent from the form of Fgeo and Fgd that the corresponding phase 
shifts A(f)geo and A0gi? do not agree with each other even if the Lorentzian term in Fgd is 
neglected. Moreover, as we will see in the next section, both phase shifts are different from 
the phase shift calculated previously in the literature as well. There are thus three different 
expressions for the phase shift for the same matter-wave interferometer. To determine the 
correct expression, we consider the low-frequency limit of these expressions, and compare 
their limiting forms to the phase shift of caused by a static Riemann curvature tensor. 

A. How the Correct Phase Shift is to be Determined 



The motion of test particles in stationary spacetimes has been well studied experimen- 
tally both directly, through measurements of the advance of the perihelion of Mercury, the 
deflection of light, and the gravitational redshift |l5| in the weak gravity limit, and indirectly, 
through the slow-down of the Taylor-Hulse pulsar 'Sli] in the strong gravity limit. It has also 
been well studied theoretically through analyses of various black hole geometries 3^. The 
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motion of test particles in nonstationary spacetimes, such as when a gravitational wave is 
present, has not been established to the same level of certainty. On an experimental level, 
until LIGO has detected a gravitational wave, there will have been no direct measurements 
of the motion of test particles in nonstationary spacetimes. On a theoretical level, while 
it would seem that the motion of test particles in nonstationary spacetimes has been well 
studied, a recent analysis [18| based on the general laboratory frame has detailed differences 
in the dynamics of test particles in nonstationary versus stationary spacetimes. 

There is, however, is no disagreement — on either on a theoretical or an experimental 
level — over the motion of text particles in stationary spacetimes. Indeed, it was explicitly 
demonstrated that in this case the general laboratory frame |la| (and thus the proper refer- 
ence frame) goes over to the standard description of the motion of test particles in stationary 
spacetimes using the geodesic EOM. 

The phase shift A(j)static of a matter-wave interferometer due to a static curvature is well 
known, and can be written as Aip static ~ niA'^T R^^^^ / h, where \A^^\ is the area of the 
interferometer, T is the transit time of the atom through the interferometer, and R^jof'^ are 
the components of the static Riemann curvature tensor in the plane of the interferometer 



[see [10], Eq. (2^6) of 



of A(l)static in 



111 ] and 



ll|, and Eq. (2.23) of 5]). While the calculation of the calculation 
5| was based on general relativity, A(f)static can be derived using only 
Newtonian gravity — where the static curvature is formed from gradients of g — as was done 
in [li|. Newtonian gravity has been experimentally verified, of course, and the effect of 
g on the phase shift of massive particles in matter- wave interferometers is well known for 
neutrons |2^, and for atoms |2^|2a]. The effect of the static Riemann curvature of the Earth 
has also been seen in measurements of the gradients in q using differential measurements of 
the phase shift caused by g at differing heights J26|, |27|. While indirect, these all serve as 
experimental indications that A(j)static holds. Indeed, there are indications that t he p hase 
shift induced by the static Riemann curvature of the Earth has already been seen (28l |. 

Consider now the limit where the period of a gravitational wave is very long compared 
to the transit time of the atom through the interferometer. The gravitational wave does not 
oscillate appreciably in the time it takes for the atom to traverse the interferometer, and to 
the atom, the curvature of the spacetime caused by the gravitational wave is effectively a 
static one. Correspondingly, the phase shift induced by the gravitational wave should have 
the same form as the phase shift induced by a static Riemann curvature. This correspondence 
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between the gravitational wave- induced, and the static-curvature-induced phase shifts serves 



M 



m 



m 



as a check of our results. When it is combined with the experimental results 

we conclude that the correct phase shift must have a low-frequency limit that has the same 

form as Acp static- 

B. The Low- frequency Limit 

Unlike Fgeo, which depends only on the transit time T of the atoms through the interfer- 
ometer, and the period 27r/u; of the gravitational wave, Fgd depends on the period Itx juj^ 
of the resonance frequency of the mirror-interferometer assembly as well. Both T and loq 
are fixed by the design of the interferometer, while lo can, in principle, take a broad range 
of values. We note, however, that the resonance frequency of the interferometer can be 
estimated as cjq ~ 2tiv sound/ L\\, where Vsound is the speed of sound of the material used to 
construct the interferometer. In almost all cases, v\\ « Vsound, and since 1/T = -^||/f||, 
27t/uo «T. 

If the transit time of the atom is much shorter than the period of the gravitational wave 
so that T << 27r/ci;, then uj << ujq as well. In this low-frequency limit, the gravitational 
wave will not oscillate in the time it takes the atom to traverse the interferometer, nor will 
the mirrors oscillate appreciably during this time. This, then, corresponds to the static-limit 
conditions considered above. 

From Eq. ^^, we see that in this low-frequency limit where uT/A << 7r/2, Fgeo ~ —3. 
Then 

iTD 771 

A(j)geo{uj) ^ —L\\L^iujhy^{u) = -3-L||L^rS^(u;), (29) 

since V^yioj) = —iujhy^{ijj)/2. Thus, /S.(t)geo{uj) is proportional to the Fourier transform of 
the Christojfel symbols for the gravitational wave in the low frequency limit, and not the 
Riemann curvature tensor. It thus has an incorrect limiting form. 

If we now take the low- frequency limit of Eq. PH| . we find that Fgd ~ —iujT/2 to lowest 
order, so that 



- 5 






A0gd(^) ~ -jrL\iL±Tu hy,{uj) = -—L\iL±RoxOy{^), (30) 

since Ro^oy^u) = uj'^h^{ijj)/2. Thus, unlike Acpgeoi^), A4>gd{^) is proportional to the Rie- 
mann curvature tensor, and has the correct limiting form. 
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The differences in the low-frequency behavior of A(f)geo and A(f)GD are shown Fig. |21 
where we have graphed A(f){Lj) / Ro^oyiuj) versus u for both expressions of the phase shift. 
For completeness, we have also included in this comparison the phase shift A(f)iTS for the 
matter-wave interferometer Fig. ^ calculated using the Linet-Tourenc-Stodolsky formula 
J3|, y]. As we mentioned in the introduction, a third calculation of the phase shift of Fig. ^ 
based also on the geodesic EOM is necessary, and in the next section this phase shift AipLTS 
is explicitly calculated for our interferometer using the Liner- Tourrenc-Stodolsky formula. 
For now, we note that the graphs in Fig. |21were plotted using the parameters of the "Earth- 
base horizontal MIGO" described pj , which has the same configuration as the interferometer 
we consider here in Fig. ^ This interferometer has a L^ = 1.2 m, L\\ = 4 km, T = 0.58 
s, Q = 1, uq = 31,000 rad/s, and used ^Li with m = 1 x lO"^'* gm. The plots in Fig. |21 
cover the same frequency range considered in [l|. Notice that for frequencies much less than 
2-71 /T =11 rad/s, A(j)GD/RoxOy approaches a constant, indicating that A(J)gd oc RoxOy for 
low frequencies, as expected, while both Acpgeo/RoxOy, and A(j)LTs/RoxOy do not. 

The behavior of all three phase shifts also differ dramatically at high frequencies. We 
will consider the high-frequency behavior of A(f)GD in appendix C. 

IV. PREVIOUS PHASE-SHIFT CALCULATIONS 

We now compare results of the above calculations with those done previously in the 
literature. 



A. Calculations done at the quantum- mechanical level 

Calculations of the phase shift by Linet and Tourrenc p, and Stodolsky ^ were both done 
at the quantum-mechanical level, but unlike us, they presented a completely general analysis, 
and did not take the nonrelativistic limit for the interfering particles they considered. Their 
expressions for the phase shift were thus manifestly covariant. Linet and Tourrenc did not 
expliatly define the coordmate system^at they used, while Stodolsky only referred to the 
presentation by Landau and Lifschitz 3.3] on synchronous reference frames. It is clear that 
since a gravitational wave in the TT gauge has h^^ = 0, the TT coordinates is a synchronous 
reference frame as well. Since both calculations were ultimately based on the action for the 
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(0 in rad/s 

FIG. 2: A plot of A(j){uj) / RoxOyi^) versus uj for A(pGD, ^(pgeo, and AcpiTS^ for the interferometer 
Fig- n using L± = 1.2 m, Ly = 4 km, T = 0.58 s, Q = 1, loq = 31,000 rad/s, and used ^Li with 
m = 1 X 10^^-^ gm. 

geodesic EOM, they were implicitly done in the TT coordinates. 

Linet and Tourrenc used the WKB approximation in the quasi-classical limit, and re- 
placed the mass-shall constraint g^^p^Pu = —mc^ (note that the signature of their metric 
is different from ours) with a differential equation in the phase of the wavefunction for the 
test particle by making the substitution p^j = {i/h)d^. Stodolsky worked directly with the 
action, and used the quasi-classical approximation we followed here to calculate the phase 
shift. The expression for the phase shift derived by both nevertheless have the same over- 
all form. The phase shift (given in Eq. [3.2.3] of l^) derived by Linet and Tourrenc does, 
however, differ from the original phase shift (given in Eq. (2.4) of |4]) derived by Stodolsky 
by a factor of 2. Since Stodolsky uses the Linet and Tourrenc phase shift expression (the 
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equation above Eq. (6.1) of [4]) to calculate the phase shift for a model interferometer, we 
will do so also. In the nonrelativistic limit of the atoms that we are interested in, this phase 
shift reduces to 

A(l)LTS = ^jj>Wh,,dt, (31) 

which is different from A0geo, even though both were calculated in TT coordinates. 

To understand the origin of the difference between ^(pLTS and Acpgeo, let us write Sgeo = 
S'gTo + Sltl where 

SlTo = \^j v'dt, Si:i = \mj v^v^K.dt, (32) 

and the superscript "kin" stands for "kinetic", while the superscript "pot" stands for "po- 
tential" . Clearly, A^g^o = A0j*^ + A^^^* . Since Mlts = Acjfgf^, contributions to the total 
phase shift A^^*" from changes to the kinetic energy of the test particle as it traverses the 
interferometer was not taken into account in Eq. (jHlj) . This kinetic term cannot be neglected, 
however. The passage of the gravitational wave through the system will cause shifts in the 
velocity of the atom as it traverses the interferometer, and thus to its kinetic energy. In 
fact, by using Eqs. (jH)) and 0, it is straightforward to show that A0^*" = — 2A0^°* + B.T., 
where B.T. are boundary terms for the atoms at the mirrors. These terms do not vanish, 
but instead contribute the 4 sine {uT /4) term to Fgeo- 

If we nonetheless use A(J)lts to calculate the phase shift for Fig. [TJ we once again find 
that A(j)LTS has the same form as Eqs. ((Hj) and (J2H1), 

In 



Flts{^T) 



sine {ujT/4:) 



2 



(33) 



The final result of our comparison of phase-shift calculations for the interferometer shown 
in Fig. [T]is given in Eq. (j^Hj) . While these equations, like Eq. (HH), are also based on the 
geodesic equation of motion using the TT coordinates. Fits nonetheless differs from Fgeo 
by the term 4 sine (u;T/4). This term comes directly from the boundary terms neglected 
in Eq. ()31|1 . (When expressed in terms of resonance functions, Eq. ()33|1 agrees, after the 
nonrelativistic limit is taken, with the phase shift calculated in j^ up to overall factors.) 

Taking the low-frequency limit, we find that F^ts ~ ~1- Thus, like Acpg^o, AcpiTS has 
the incorrect low-frequency limiting form. It also has a /iz^^/i-frequency behavior that is 

22 



dramatically different from either Fgeo or Fgd-, as can be seen in Fig. |21 In this limit, while 
Fgeo ~ \/ujT and Fgd ~ ^J.^^^ appendix C), Flts ~ l/(^^)^5 ^^^ dies off much more 
rapidly. As pointed out in [l|, this rapid die off is due to three effects. Firstly, the A(f)iTS 
does not include the boundary terms at the mirror. These terms depend on the value of the 
velocity of the atom at the surface of the mirror, and the length L\\/2, and they mitigate 
this rapid decrease for A0geo- Secondly, AcpLTS depends on v^, and is constant as the atom 
travels through the interferometer, while A(J)gd depends on the product xv, and increases 
linearly with the distance the atom travels from the initial beam splitter. Thirdly, since 
AcpLTS depends only on v'^v^hij, discontinuous jumps in the velocity due to the hard wall 
mirrors were not taken into account in A(f)LTS- From appendix C, these hard-wall jumps 
contribute directly to the high-frequency behavior of A(J)gd- 

B. Calculations done at the quantum-field-theory level 

1. For Spin Fields 

Calculations of the phase shift at the quantum-field-therapy level were first done by 
Cai and Papini jg]. Their starting point is the covariant version of the Ginzburg-Landau 
equation, but for a charged particle. This is a relativistic Klein-Gordon equation for a 
bosonic field Lp that interacts with itself through a Iv?]^ interaction (Eq. (2.1) of [q|). By 
turning off the \ip\'^ interaction, and setting the charge of the particle to zero, one recovers 
Linet and Tourrenc's starting point. If one extracts the single-particle Hamiltonian from 
this equation, and take the correspondence principle limit of the EOM, one recovers the 
geodesic EOM used by Stodolsky. Thus, like Linet and Tourrenc, and Stodolsky, Cai and 
Papini work in TT coordinates. 

By acting on if with a phase operator x, Cai and Papini was able to solve for a general 
expression for the phase shift due to gravitational effects in the weak-gravity limit through 
a stationary-phase type of analysis. This operator has the form [Eq. (3.1) of |j5|]], 

X=[-\jdz^ (r«,A/3 - r^,,„) r^ + \j dz^g^pPf\ , (34) 

(using our notation) for a closed path where Jap is the generator of Lorentz transformations, 
and P" is the generator of translations. We have dropped the vector potential term since the 



23 



atoms we consider in this paper are not charged. As they mentioned, using Stokes' theorem 
the first term in x can be written as 

-7 / R,uo.pr^dT^\ (35) 

4 Jv 
in the weak-gravity hmit; the integral is over the timehne surface V bounded by 7, and 
dr^^ is a 2-form on V. This expression has been exphcitly verified in lUJ] for the phase shift 
caused by Riemann curvature tensor for a stationary spacetime. 

The second term agrees with Stodolsky's expression for the phase shift. Indeed, Cai and 
Papini uses Eq. flMj) to calculate the phase shift of a model interferometer due to the passage 
of a gravitational wave [Eq. (3.15) of [5|], and they found that up to terms of order (f /c)^, 
the following holds: 

"The Ax2 terms in (3.15) agrees with Stodolsky's result [11] when a misprint in 
the latter is corrected by inserting an overall factor |." (From page 415 of J5|-) 

The reference [11] cited above is the work by Stodolsky cited as J| in this paper. Con- 
sequently, their expression for the phase shift, like Linet and Tourrenc's, and Stodolsky's, 
also has the incorrect low-frequency limiting form. The effects of the mirrors on the test 
particle's phase was not considered by them either. 

2. For Spin 1/2 Fields 



Borde and coworkers ja, LD, |8| considered the use of fermionic atoms in matter- wave inter- 
ferometers to measure various gravitational effects, including the detection of gravitational 
waves. Their starting point was the Dirac equation for a two- level, spin 1/2 atom expressed 
in curved spacetimes using the tetrad formalism. Similar to Cai and Papini, theirs was a 
Hamiltonian-based approach, and they showed that the dominant contribution by the grav- 
itational wave to the Hamiltonian is from a term with the form piPjh^^ /2m, where m is 
the mass of the atom, and pi its momentum [see Eq. (10) of 6]]. The other terms in the 
Hamiltonian due to the gravitational wave are from very small spin-gravitation couplings. 
These terms can be neglected here since we are dealing with atoms without spin. When 
this is done, the Hamiltonian they calculated has the same form as Eq. (J2I) for the geodesic 
EOM, and is manifestly different from the Hamiltonian Eq. p7|) for the geodesic deviation 
EOM. Thus they, too, worked in the TT coordinates. 
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Although an expression for the phase shift was not calculated in p, they concluded the 
following: 

"The phase shift, that one can calculate from the first terms of Eq. (10), increases 
with the velocity of the atoms in contrast with what happens in the inertial 
field case and a gravitational wave detector using atomic interferometry should 
therefore use relativistic atoms. In this case, it can be shown directly from Eq. (5) 
that the phase shift is proportional to /i(x,+) / c?s/A where A is the atomic de 
Broglie wavelength, result which is similar to the optical case." (From page 161 



of 



!•) 



In a subsequent paper J8|, Borde gave an explicit expression for the phase shift [see Eq. (98) 
of J8|]. It has, for its spin component, the same expression, A^^y^, for its phase shift as 
Stodolsky's. Thus, like Stodolsky's approach, any phase shift calculated using this analysis 
would have the incorrect low-frequency limiting form as well. 

3. More General Analysis 

Finally, all of the analyses that we outlined in this paper were done in the weak-gravity 
limit. In J2|, a derivation of the phase shift of a quantum mechanical test particle was done 
without making this approximation. Their analysis was also based on the geodesic-EOM 
approach, and they arrived same expression for the phase shift that Stodolsky does for spin 
1/2 test particle with no approximations, while it holds only the first order in h for spin 
and spin 1 particles. Consequently, phase shifts calculated with their expression will also 
have the same incorrect low-frequency limiting form as Stodolsky's. 

V. CONCLUDING REMARKS 

There has been a preference among those who study matter-wave interferometry and 
gravitational waves to use TT coordinates to analyze the phase shift of an interferometer 
induced by gravitational waves instead of the proper reference frame. The expectation is that 
irrespective of which reference frame is used to calculate the phase shift, the final expression 



will be the same 
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23[. We have shown in this paper that for matter- wave interferometry 



this expectation is wrong. By presenting a calculation of the phase shift using the TT 
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coordinates following the line of the analysis found in the literature 3|, y, |5|, ly, lD, |8|, |9| , and 

also a calculation of the phase shift using the proper reference frame following the analysis 
found in l|,y], we have shown that the expressions for the phase shift obtained using the TT- 
coordinates is different from that obtained using the proper reference frame. Moreover, we 



have shown that the commonly accepted expression for the phase shift Eq. 



which is the 

is 






end result of all previous derivations of the phase shift in the literature 

incomplete, and gives yet a different expression for the phase shift of Fig. ^ This happens 

even though derivations of this expression were all done in TT coordinates. 

The focus of this paper has been on demonstrating in detail that differences exist between 
the various approaches to calculating the phase shift, and showing explicitly that different 
expressions for the phase shift are obtained. To determine which one, if any, is correct, we 
noted that in the limit of sufficiently low frequency gravitational waves, the correct expression 
for the phase shift must reduce to the well known expression for the phase shift induced by a 
static curvature. It must therefore be proportional to components of the Riemann curvature 
tensor. This static phase shift can also be calculated within Newtonian gravity, where the 
curvature is formed from gradients of the acceleration due to gravity g. The phase shift due 



to g has been measured yA 



well 
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25l |. and recently has been used to measure gradients in g as 



28[. In addition, there are indications that the phase shift induced by the static 
Riemann curvature tensor has already been seen [28|. We have explicitly shown that when 
the phase shift is calculated using the proper reference frame and geodesic deviation EOM, 
it has the expected limiting form in the low- frequency limit, and is therefore correct. The 
other two, both of which used the TT coordinates and the geodesic EOM, do not, and are 
thus incorrect. 

The implications of this conclusion are immense. In [l|, y], MIGO, the Matter- wave 
/nterferometric Gravitational-wave Observatory, was proposed. The theoretical sensitivity 
predicted for MIGO was based on a calculation of its phase shift using the proper reference 
frame and the geodesic deviation EOM. The conclusion arrived at in these papers was that 
by using slow atoms, it would be possible to construct gravitational wave observatories that 
are as sensitive as LIGO or LISA (the Laser Interferometer S'pace Antenna), but would be 
orders of magnitude smaller, and have a much broader frequency response. This conclusion 
is diametrically opposite from the conclusions arrived previous in the literature that fast 
atoms — indeed, ultrarelativistic ones — should be used in the matter-wave interferometer 
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instead p. 

The conclusion that relativistic atoms should be used in the detection of gravitational 
waves with matter-wave interferometers is the basis on which most of the atomic, molecular 
and optical physics community judges the suitability of using matter- wave interferometers 
as gravitational wave detectors. The physical origin of this conclusion by Borde and his 
coworkers is the dependence of the phase shift AipiTS on the square of the atom's velocity 
in the geodesic-EOM-based approach; clearly, the faster the atom, the larger the phase shift 
will be. Even though our analysis of the phase shift was done for nonrelativistic atoms, 
we note that if the size of the interferometer is held constant, and only the velocity of the 
atom is increased, then in the high-velocity limit, T ^ 0, and eventually ujT « 7r/4 for the 
frequency of gravitational waves expected from astrophysical sources. The high- velocity limit 
is thus equivalent to the low-frequency limit for the gravitational waves, and we have shown 
that the low-frequency limit of A(f)LTS is incorrect. Thus, Borde and coworker's assertion 
about relativistic atoms being more suitable for the detection of gravitational waves is also 
incorrect. 

All previous calculations and estimates of the phase shift in the literature also neglected 
the effect the mirrors of the interferometer have on the phase shift. Indeed, from appendix 
C we see that it is precisely because hard wall mirrors are used — along with the instanta- 
neous impulsive force that they exert on the atoms when they reflect off — that the use of 
slow-velocity atoms leads to a sensitivity that increases with frequ ency is obtained. This 
dependence on frequency is the exact opposite for LIGO and LISA ^], and is part of the 
reason that MIGO can be so much smaller than either one. 

We have not, in this paper, presented the underlying physical reasons why the current 
TT coordinates approach to calculating the phase shift fails, and the proper reference frame 
succeeds. We shall address this issue in detail later J3J|. Here, we only note that one may 
try to modify the boundary conditions for the atom on the mirrors so that the phase shift 
calculated for Fig. ^in TT coordinates will have the correct low-frequency limit. Doing this 
is especially tempting in light of the following statement by Thorne: 

"However, errors are much more likely in the TT analysis than in the proper- 
reference-frame analysis, because our physical intuition about how experimental 
apparatus behaves is proper-reference-frame based rather than TT-coordinate 
based. As an example, we intuitively assume that if a microwave cavity is 
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rigid, its walls will reside at fixed coordinate locations x^ . This remains true 
in the detector's proper reference frame (aside from fractional changes of or- 
der {L^/X )h, which are truly negligible if the detector is small and which the 
proper-reference- frame analysis ignores). But it is not true in TT coordinates; 
there the coordinate locations of a rigid wall are disturbed by fractional amounts 
of order h, which are crucial in analyses of microwave-cavity-based gravity wave 
detectors." (From page 400 of |l2l|.) 

We find that while it may be possible to suitably modify the boundary conditions to bring 
A^peo in line with the low frequency limit, doing so will be ad hoc, and will not explain the 
underlying physical reason for the differences between A(j)geo and A(J)gd- What is at the 
heart of the physics are not questions of which boundary conditions to take and how they 
should be modified, but rather how measurement of physical quantities are made in general 
relativity when the background metric varies with time. 

When doing an experiment to measure some physical property of a test particle — its 
phase shift, say — we start by choosing an origin for our coordinate system. It is natural 
to choose as this origin some point on our measuring apparatus, which in our case is the 
initial beam splitter. In the three spatial dimensions, this origin is represented by a point 
on a three-dimension space; this point is what we use when measurements with a physical 
apparatus are actually made. If we wish to represent the origin in space and time, it would 
seem to be natural to represent it as an event in four- dimensional spacetime. Coordinate 
systems constructed by doing so are called quasi-cartesian coordinate systems by Synge. He 
also noted who also noted that a suitable choice of origin in spacetime for is often difficult 
(see section 10 of [llf), however. This is especially true in the case of gravitational waves, 
where the metric is nonstationary, and varies with time. In fact, by identifying the origin 
as an event, we have ignored the fact that our measuring apparatus is a physical object. It 
too must travel along a worldline in the spacetime, as will any test particle whose physical 



properties we are measuring. As was argued in [18[ , only relative motions of a test particle 
with respect to the apparatus can thus be measured. These arguments were based in part 
on the following observation by Hawking and Ellis: 

"In chapter 3 we saw that if the metric was static there was a relation between the 
magnitude of the timelike Killing vector and the Newtonian potential. One was 
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able to tell whether a body was in a gravitational field by whether, if released 
from rest, it would accelerate with respect to the static frame defined by the 
Killing vector. However, in general, space-time will not have any Killing vectors. 
Thus one will not have any special frame against which to measure acceleration; 
the best one can do is to take two bodies close together and measure their relative 



acceleration [emphasis added]." (From page 78 of 



m- 



This relative motion is not measured in the TT coordinates where the coordinates are defined 
relative to a fixed event, chosen as the origin, in the spacetime. The motion of the measuring 
apparatus through spacetime was not taking into account in the TT coordinates, so that the 
beam splitters and the mirrors are at fixed coordinate values. This relative motion, on the 
other hand, is measured in the proper reference frame where the coordinates are explicitly 
constructed to span the distance between the measuring apparatus, and the test particle 
whose properties are to be measured. The motion of the measuring apparatus has been 
built in from the beginning, and indeed, the positions of the beam splitters and the mirrors 
shift during the passage of the gravitational wave. Failing to take into account the motion 
of the measuring apparatus in spacetime is the underlying physical reason why the wrong 
expression was obtain for Fig. ^ when the TT coordinates are used. 
These statements will be expanded upon, and justified in |34] . 

APPENDIX A: THE QUASI-CLASSICAL APPROXIMATION AND THE IN- 
TERFEROMETER 

In this appendix we set the framework used in our different calculations of the phase 
shift. 

1. The quasi-classical approximation 

Like Stodolsky y, |9|, we use the quasi-classical approximation to calculate the phase 
shift of an atom in a matter-wave interferometer. This approximation is the Lagrangian 
version of the WKB, or stationary phase approximation, in the Schrodinger representation 
of quantum mechanics. We assume that the density of particles through the interferometer 
is low enough that they do not interact with one another appreciably in the time it takes 
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for the atom to traverse the interferometer. Within this approximation, the particles can 
be treated as test particles that interact only with the gravitational wave, and not with 
themselves. We can thus consider the phase shift of each atom individually, and use the 
Feynman path integral method to calculate its phase shift. However, instead of integrating 
over all possible paths, we make the quasi-classical approximation, and consider only the 
most probable path — which is also the classical path — that the test particle takes through 
the interferometer (see |36j). As a result, 

A0^ ^lSci[xA,XB;-fi]- S^^[xa,xb;12] \, (Al) 

where the superscript "cl" stands for "classical" , and 71 and 72 are the two possible paths a 
classical test particle can take through the interferometer that links a common point xa on 
the interferometer (the initial beam splitter, say) to a common point xb on the interferometer 
(the final beam splitter, say). Note that Sci[xA,XB]'y] are functionals of the points xa and 
Xb, along with the paths 71 and 72 connecting them. 

Unlike ,4], as well as the calculations of the phase shift done in |3|, |5l ly, [D, |8|, |9| , we will 
work in the nonrelativistic limit. The reason for doing so is two fold. First, the proper 
reference frame, as well as the general laboratory frame and all the other reference frames 
based on the worldline of an observer, are constructed for nonrelativistic test particles. 
Second, as is well known, the phase shift of a matter-wave interferometer is proportional to 
the mass of the particle. This argues for the use of atoms in matter-wave interferometers 
to detect gravitational waves, as was proposed in jl|, |2] , and in any conceivable matter- wave 
interferometer, these atoms will be moving nonrelativistically. 

2. The interferometer 

The configuration of the interferometer that we use in our calculations is shown in Fig. d 
Details of this interferometer can be found in [l|. For our purposes, it is sufficient to note 
that the interfering test particle in the interferometer is an atom with mass m. A beam 
of these atoms is emitted from an atomic source with velocity Vg, and is beam-split into 
different diffraction orders using a diffraction grating with periodicity a (see J38l |39l |40|). 
Only the n = +1 and n = — 1 diffraction orders are used in the interferometer. 

There are two possible paths, denoted as n = +1 and n = —lin Fig.^ that any individual 
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atom can follow through the interferometer; it is not possible to tell a priori which of the 
two paths an atom will take. If the atom is diffracted into the n = +1 diffraction order, then 
immediately after the beam splitter it will have a velocity of +v_[_ = 27rh/ma perpendicular 
to the longitudinal axis of the interferometer. Similarly, if an atom is diffracted into the 
n = —1 diffraction order, it will have a velocity of —v± (see Fig. Q). Since the diffraction 
grating does not impart energy to the atom, for both diffraction orders v\\ = (ff — v'^Y^'^ is 
the velocity of the atom parallel to the axis of the interferometer immediately after it leaves 
the beam splitter. 

In n a .he use of haM-waai. specula. .eflec.o„ „..o. w. con.ae.ed f„. .he .a.e. 
interferometer. Mirrors of this type have been made from silicon wafers [33] , and for near 
normal and near glancing angles, the reflection from these mirrors is close to 100%. These 
are not the mirrors that are currently used in matter- wave interferometry, however. Mir- 
rors based on diffraction gra ting s, standing waves of light, or stimulated Raman pulses are 
271 . l38l . I39L l40l | . We will nevertheless consider only the case of hard- wall 
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mirrors in this paper, because of their simplicity. In the absence of gravitational waves, the 
effect of hard-wall mirrors, if they do not move, on the velocity of the atoms through the 
interferometer is 

Vx{out) = Vx{in), Vy{out) = —Vy{in), (A2) 

where v{in) and v{out) are, respectively, the velocity of the atom immediately before it is 
reflected off the mirror, and immediately afterwards. 

APPENDIX B: THE TT GAUGE AND THE GENERAL LABORATORY 
FRAME 

In this appendix we will show how the TT gauge can be taken in the general laboratory 
frame in a treatment that follows (l8l |. 

Like the proper reference frame, the general laboratory frame is an orthonormal coordi- 
nate system fixed on the worldline of an observer. In our case it is fixed on the worldline 
of the initial beam splitter in Fig. ^ We choose a local tetrad e^ fixed on the observer's 
worldline. Thus, e"^ ^cau = g^iu, while e^'^esu = Vab- In a slight break in the notation 
followed elsewhere in this paper, which follows the notational convention established in lial , 

n 

we use in this appendix the convention taken in [18||. Thus, we use as the tetrad indices 
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capital Latin letters, which run from to 3. Lower-case Latin letters are used as the spatial 
tetrad indices, and they run from 1 to 3. As usual, Greek letters are used for the spacetime 
indices of a general coordinate x'^. 

Consider a the motion of a test particle V close to the observer O and her worldline. Let 
X^{t) be the position of V at any proper time r in the observer's frame. X^{x) : x^ -^ X^ 
can also be considered as a coordinate transformation from the general coordinates x^ to 
the tetrad frame at any time r. Thus in a small neighborhood Uq of O, 

dX^ 



dX^ 
r]AB dX^dX^ = r]AB-Q^ 

so that 



Uo ^^" 



dX' 



A dx^ 



dx^^ 
Integrating Eq. ()B2|) . we find for the spatial coordinates 



dx^'dx'', (Bl) 



(B2) 

Uo 



Xa{r)= /' ' e„^(r,s)x^(r,s)c/s= / e„. (B3) 

where e^ = CAfj.dx'^ is a 1-form, and 7^ is a space-like geodesic linking O to V at a given 
proper time r of the observer. Next, after taking appropriate account of causality, the time 
component is 



Xo(t) = -(r-r')+ / eo. (B4) 

J 7^(0-) 

where 7,^ is a null-vector linking O at time r' with V at time r; the difference r — r' is the 
transit time of a light pulse emitted by O, and scattered by V. 

In the case of linearized gravity, g^^i, = rj^y + h^y. It is straightforward to see that in terms 
of tetrads 



K-,,K]8a. and e^,= {K + ,:K]5t. (B5) 



^/^e)^A, and e\=[5; + \..^ 

This choice is not unique. If one does a local Lorentz transformation L^"^ on e^'^, then 
TjAB = L^^L^^ while e^'^ = L^'^e^'^. This leaves g^i, = ca^i ^^ y = ^au ^^ v invariant. 

The usual derivation of the TT gauge makes use of a remnant of general coordinate 
transformation invariance still left in the linearized-gravity limit of general relativity. If 
2^M ~ ^M + ^M where ^^ is a 'small', arbitrary vector, then h^^y transforms to 

~hi,v = h^y + d^iy + dyii,. (B6) 
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One can choose ^^ such that 

can be chosen for GWs in Minkowski space. Doing so defines a specific coordinate system, 



which is called the TT coordinates in 



12 



The TT gauge can also be implemented in the general laboratory frame as well. In terms 
of a tetrad frame, the linearized coordinate transformation Eq. ()B6|1 induces a Lorentz 
transformation on the tetrads such that e"^ ^ = L^^c^^, where 

L^A = ^^A + l {d^U + Oa^^) , (B8) 

and Oa = e^c?^. Like the TT coordinates, the TT gauge can also be realized, but now 
through the appropriate Lorentz transformation. From Eq. ()B2|) . this in turn means we can 
choose the TT gauge in the construction of the general laboratory frame as well. 



APPENDIX C: THE HIGH-FREQUENCY LIMIT 

In this appendix we consider the high-frequency limit of A0gd- 

Because of its dependence on three timescales, the analysis of the high-frequency limit 
of A(f)GD is niore involved than that for the low frequency limit. There are two cases to 
consider. The simplest case is when the transit time of the atom through the interferometer 
is very long in comparison with the period of the gravitational wave, so that T » 27t/uj, 
but u << ujq still. To the atom, the mirrors still do not move in response to the gravitational 
wave, but the atom will experience many oscillations of the gravitational wave as it traverses 
the interferometer. In this limit the phase of the atom oscillates so fast that we would naively 
expect the phase shift to average to zero. It does not do so, however. We instead find that 

Fgd ~ 1- 

That A(j)GD does not vanish in the limit of high-frequency gravitational waves is due to 
action of the mirrors on the atoms as they traverse through the interferometer. When the 



atoms refiect o 
instantaneous 



le mirrors, they impart a impulsive force on the atoms that is effectively 
4l|. The effect of this force is readily apparent from the discontinuity in Vy 
across the mirror in Eqs. ()A2|) and ((Tj). Note in particular that the force is proportional 
to perturbed velocity viy{tr + T^^/2) of the atom caused by the gravitational wave, and is 
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different depending on which path the atom travels through the interferometer, as can be 
seen from the force hnes in Fig. ^ It is because of this net instantaneous force that A0gd 
does not average to zero in the high-frequency hmit. 

When, uj » ujq as well, the mirrors themselves oscillate rapidly as the atom traverses 
the interferometer. This rapid motion of the mirror is not correlated with the motion of the 
atoms through the interferometer, however. When the atom reaches the mirror, roughly half 
the time it will be moving in the same direction as the mirror, which would, from Eq. (fTHjl . 
tend to decelerate the atom. The other half of the time it will be moving in the opposite 
direction as the mirror, and would tend to accelerate the atom. Not surprising, in this limit 
the resonance function decreases by a half, and Fgd ~ 1/2. 

In this high-frequency limit A(J)gd is T^-ot proportional to RoxOy There is no reason to 
expect it to be, as there was in the low-frequency limit, however. Because the mirrors exert 
a net force — which itself is the result of the passage of the gravitational wave through the 
system — on the atom as it passes through the interferometer, Acpcn ~ Toy 
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